Universality in nonequilibrium condensation of exciton-polaritons 
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We investigate the process of condensation of exciton-polaritons in a one-dimensional nanowire, 
predicting spontaneous formation of domains of uncondensed excitons and condensed polaritons. We 
demonstrate a universal scaling law for the density of domains, which results from the competition 
between characteristic timescales present in the system. However, we find that the system does not 
follow the standard Kibble- Zurek scenario of a nonequilibrium phase transition. 
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Semiconductor microcavities, in which photons are 
confined and strongly coupled to excitons, allow for cre- 
ation of exciton-polaritons. These bosonic quasi-particles 
can undergo Bose-Einstein condensation at standard 
cryogenic or even room temperature [H-Q. The achieve- 
ment of exciton-polariton condensate opened a path for 
new technological developements, eg. in optoelectronics 
(ultra fast switches, quantum circuits) [5[ or medicine 
(compact terahertz lasers) Q. It is known that condensa- 
tion can have a character of a nonequilibrium phase tran- 
sition when it takes place on a finite timescale Q. Such 
transitions have been attracting great attention in many 
fields of science Q. One of the early studies concerned 
topological defect formation in cosmological models [jj. 
As the hot universe cooled down, the phase transitions 
broke symmetries of the vacuum fields, which leaded to 
the creation of topological defects. This mechanism has 
been also considered in a wide variety of other physical 
models, including superfluid Helium (To| . superconduc- 
tors, cold atomic gases, and other systems (llj. The out- 
come of the transition is the creation of defects such as 
domain walls, monopoles, strings, vortices or solitons 0- 
Bose-Einstein condensates of exciton polaritons, as 
highly controllable nonequilibrium quantum systems, are 
very suitable for studying these phenomena. 



The Kibble-Zurek mechanism (KZM) theory is a pow- 
erful tool that allows to predict the density of defects 
without solving the full dynamical equations. It is as- 
sumed that the system adiabatically follows the equilib- 
rium state up to the "freeze-out" time t before the crit- 
ical point. At this moment the relaxation time becomes 
too long for the system to adjust to the change of pa- 
rameters. Consequently, the fluctuations approximately 
freeze, and the system enters the "impulse" phase. Af- 
ter crossing the critical point, the adiabaticity is restored 
when the relaxation time becomes short enough again. If 
the system in the final phase exhibits symmetry break- 
ing, distant parts will choose to break the symmetry in 
different ways. The number of resulting defects will be 
determined by the "memorized" healing length £ at the 
time t. Remarkably, the time t and the number of defects 
iVd can be determined solely from the critical exponents 



of the phase transition, according to the scaling laws 

|£| ~ „-*"/(!+*") ; iv d - £~ d ~ V du/{1+Zl,} , (1) 

where v is the rate of parameter change, d is the dimen- 
sionality, z and v are the critical exponents. All the de- 
tails of complicated dynamics are hidden, but we are still 
able to derive Eqs. ([TJ thanks to the underlying universal- 
ity. The Kibble-Zurek scaling laws have been confirmed 
numerically in many cases, including the condensation of 
ultracold atoms 0. 

The condensation of exciton-polaritons, on the other 
hand, is an example of a qualitatively different case. This 
phase transition may be named "nonequilibrium" in yet 
another meaning of this word. In fact, there may not 
be a well defined equilibrium state that the system could 
follow, due to the lossy nature of the system itself 0. 
The particle loss and gain may be so strong that they to 
a large extent determine the relaxation of the system. 

A natural question arises, whether the Kibble-Zurek 
scenario does apply to the case of an intrinsically 
nonequilibrium system such as the polariton condensate? 
And do the critical scaling laws such as the one given by 
Eq. ([1]) hold? In this letter, we aim to answer these ques- 
tions by investigating the exciton-polariton domain for- 
mation during condensation in a one-dimensional GaAs 
nanowire We note that recently it was suggested 
that the KZM describes the formation of quantum vor- 
tices in a related two-dimensional system [J]. However, 
no evidence of the above adiabatic-impulse-adiabatic sce- 
nario or the critical scalings was given. We show that in 
the case considered here the system does not follow this 
scenario, and consequently the dynamics cannot be ad- 
equately described by the standard KZM. Nevertheless, 
power laws similar to Eq. ([I]), originating from the com- 
petition of characteristic timescales in the system, are 
still present. 

We model polariton condensation in an incoherently 
pumped nanowire within the truncated Wigner approxi- 
mation adapted to the polariton case ■ The stochastic 
equations for the polariton order parameter tp(x,t) and 
the reservoir exciton density tir.(x, t) are (in the quasi- 
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one-dimensional version) 



dip 



-dt 



ip + dW, 



(2) 



dn R 
dt 



= P(x,t) - lR n R - R w n R \^\ 



where m* is the effective mass of lower polaritons, U (x) is 
the disorder potential due to imperfections of the sample, 
R is the rate of stimulated scattering, 7c and j R are 
the polariton and exciton loss rates, gc and g R are the 
respective interaction coefficients, P(x, t) is the external 
pump profile, and dW is a complex stochastic variable 
where 



(dW{x)dW{x')) =0, 

(dW(x)dW*(x')) = ^-(R 1D n R + lc )5 x , x , , 
2Ax 



(3) 



reflecting the quantum noise due to the particles entering 
and leaving the condensate. The coefficients R, gc, g R 
are renormalized in the quasi-one-dimensional case by 
assuming a Gaussian perpendicular profile of \ip\ 2 and 
n R of width d determined by the nanowire thickness, 
(i? 1D , 5^ 1D ) = ( J R, ff^ )/^/2^. 

We consider the simplest case of a flat pump beam 
of extent L, switched on at the time t = 0, P(x,t) = 
Po6(t)9(L/2 — \x\), with 8(x) being Heaviside step func- 
tion. Nevertheless, the general conclusions should hold 
for other spatial and temporal pulse shapes, as well as 
electrical pumping, as long as the pump is switched on 
rapidly. The above model describes the dynamics of the 
system at zero temperature. In reality, the distribution 
of excited modes may resemble a thermal state with tem- 
perature comparable to the ambient temperature 0, [l3| . 
However, such mode distribution should not influence the 
phenomena presented here. We will see that the corre- 
lation length in the condensed phase is equal to several 
tens of micrometers, order of magnitude larger than the 
correlation (de Broglie) length in the thermal state at a 
typical temperature of 10 K. Therefore, the excited mode 
distribution can be well approximated by a flat distribu- 
tion in momentum space, in accordance with Eqs. ([3]). 

In Fig. [1] we present results of a typical simulation of 
the condensation process. The FigurcsQJa) and (b) show 
the polariton and reservoir exciton density in a single sim- 
ulation, which can be interpreted as a single experimental 
realization (pulse). The polariton condensate appears af- 
ter about 350 ps after switching the pump on, in the form 
of areas of high density, complemented by areas of high 
incoherent exciton density. While the polariton domains 
appear to be vulnerable and erratic structures (dissolving 
after few hundreds of picoseconds), they become much 
more pronounced and well defined when averaged over 
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FIG. 1. Evolution of (a) polariton density \tp\ 2 and (b) exciton 
reservoir density nn during condensation in a ID wire with 
a disorder potential of amplitude Uo = 0.05 meV using a top 
hat beam with Po = 5 Pth. (c) Polariton density averaged over 
100 pulses. The magnified image shows three of the visible 
polariton domains created due to the nonequilibrium nature 
of the phase transition. Other parameters are 7c = 1/30 ps, 
7R = 1/400 ps, d= 5^m, R = 2.8xlO"Vm 2 ps" 1 , g c = SfJ-eV 
fim 2 , g n = 2^0- 



many realizations of the stochastic equations, which cor- 
responds to averaging over many pulses, see Fig. QJc). 
Here we included a small realistic disorder potential U(x) 
with amplitude of fraction of meV and correlation length 
of a few micrometers. It allows to pin the domains in 
certain regions of the sample, that would otherwise be- 
come smeared out over many realizations. The apparent 
exciton-polariton phase separation originates from the re- 
pulsion between the two components. It is related to the 
existence of two stationary solutions to the homogeneous 
version of Eqs. ©. Substituting ip{x,t) = ipoe'^ and 
n R (x,t) = n R0 we get 



(a) |^ | 2 =0, 

(b) Ho\ 2 =V 



7R 



P 1D ' 



"-R0 = — , (4) 
7R 
7c 

HR0 - -^Jd , 



where ij = Po/Pth — 1 and P t h = 7c7r/-R 1D is the thresh- 
old value of the pump power. These are the exciton and 
polariton phases, the latter existing only for Pq > P t h- 
Above Pth the system may choose locally to approach 
one of them. 

Let us turn our attention to the initial dynamics of the 
system. Initially, both the polariton and exciton modes 
are almost empty. After switching the pump on, the 
reservoir exciton density n R in the illuminated area starts 
to increase gradually. The polariton field, on the other 
hand, is damped or amplified depending on the value of 
e(x, t) = R m nR — 7c, see Eq. ([2|). The condition e = 
determines the critical point, and e has the same role as 
the relative temperature in the case of a thermal phase 
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FIG. 2. Evolution of the mean polariton density \i/j\ 2 and 
reservoir density tir in Fig. [T] averaged over the illuminated 
area. The vertical line corresponds to the condensation 
threshold of iir when the polariton component becomes un- 
stable. The dotted lines correspond to the analytical re- 
sults ©-J6]). 



transition [lfj . We are not able to control e directly, but 
by varying the pump power Pq, we can control the rate 
of the phase transition by changing the slope of e{i) at 
the critical point. 

The initial evolution of polariton and exciton occupa- 
tions can be well approximated with an analytical solu- 
tion which neglects the nonlinear terms in Eqs. ((2]). For 
convenience, we now choose t = at the time where the 
reservoir exciton density reaches the polariton amplifica- 
tion threshold jir = n t h = R 1B /jc, and obtain 



7c 
R 



Pth 



■qe 



-7R* 



N(t) = N(0) exp 



7c 
7R 



V (iRt 



-7R* 



1) 



(5) 



(6) 



where N(t) = J \ip(x,t)\ 2 dx is the number of condensed 
polaritons. Figure [5] shows excellent agreement of this 
simplified model with numerics up to the instant when 
the polariton density saturates due to the nonlinearity. 

The crucial observation is that while polariton density 
\ip\ 2 remains small, the relaxation is completely deter- 
mined by the gain and damping terms, and the related 
timescale is given by f re ] = l/\e\. One of the conse- 
quences of strong noncquilibrium nature of the system is 
that there is practically no mixing of modes with differ- 
ent momentum, and consequently the relaxation does not 
establish any relevant length scale (healing length), as 
e = 7c?7 (1 — e m ) is independent of momentum. Con- 
trary to the KZM, there is no well defined equilibrium 
before the phase transition, and the mode occupation is 
mainly shaped by the noise correlations ([3]). Moreover, 
as we show below, it is not these initial correlations that 
determine the final number of created defects. 

The situation changes when the polariton-polariton in- 
teraction becomes sufficiently strong to compete with 
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FIG. 3. Average number of polariton domains as a function 
of the relative power in bilogarithmic scale. Full points corre- 
spond to experiments averaged over 30 different samples with 
Uo = 0.05 meV, and empty points correspond to a sample 
without disorder. For each sample, 50 Wigner realizations 
were done. The solid line is the solution of the analytical 
model Eqs. l[5 | -([8 | . The scaling exponents from linear fits 
(dashed lines) in the case 7c = 1/30 ps are 0.24 (with disor- 
der) and 0.29 (no disorder), and the respective exponents in 
the case 7c = 1/120 ps are 0.25 and 0.36. 



noncquilibrium processes 

e(i) = ag^n(i) , 



(7) 



with a being a constant of the order of unity, and 
n(t) = N{t)/L the average polariton density in the il- 
luminated area. At this point t, which can be identi- 
fied as the condensation time [13j, mode mixing starts 
to take place, which results in the appearance of an in- 
stantaneous quasi-equilibrium. It is this moment, when 
the length scale associated with the size of emerging do- 
mains is established. The number of created domains is 
determined by the interaction-dependent healing length 
e= (h 2 /2m*g}Pny 



1 /2 

m' g^Pn) at the time £, 



N d ~L/S~(g%> n(t)) 



a/2 



(8) 



with L being the system size. 

In Figure [3] we show the dependence of the average 
number of created domains [Tij on the pumping parame- 
ter r]. In the absence of disorder, the agreement between 
the numerical data and the theoretical model of Eqs. (0)- 
([8]), depicted by solid lines, is very good over the power 
range of Po = 2P t h ■ • ■ 20Pth for two different values of 
7c- This confirms the validity of our scenario of do- 
main formation, based on the transition from the strongly 
nonequilibrium to the quasi-equilibrium regime. We note 
that the depicted dependence is well approximated by 
power-law functions (straight dashed lines) both in the 
presence and absence of the small disorder potential. In 
the presence of disorder, the slope is somewhat smaller, 
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FIG. 4. Theoretical predictions for the number of domains 
Nd as a function of rj. Solid line corresponds to the analyt- 
ical approximation ([9]), dotted line is the numerical solution 
of Eqs. ©-([H]), green dashed line corresponds to the limit 
V 3> Tr/yc, and blue dashed line to the limit rj <C Yr/7c- 
Parameters are the same as in Fig. with 7r = 1/30 ps. 

which we contribute to the ability of disorder wells and 
hills to shape additional domains. We found out that 
even in the case of very small Pq, there was a minimum 
number of domains induced by the presence of the disor- 
der. 

A simplified analytical solution can be obtained by ex- 
panding Eq. ([5]) up to the second order in jut. Doing 
this, one deals with third-order equation for the conden- 
sation time t which has the physical (real and positive) 
solution 

2YRi = l + ^ + £, (9) 

where D = (a-1+ ^Ja(a - 2)) 1 / 3 and a = I2C7H./777C 
with C logarithmically dependent on the parameters. 
The resulting scaling of domain number is presented in 
Fig. SI In the limit rj 3> 7r/yc the number of domains 
scales like Nd ~ rj 1 / 4 , while in the limit rj <C Yr/yc the 
scaling is N4 ~ 77 1 / 2 . The transition between the two 
regimes is determined roughly by 77 ss yr/YCi an d takes 
place over few orders of magnitude of rj. Therefore, for 
a limited range of 77, the scaling law may still appear as 
power- law with exponent between 1/2 and 1/4, in agree- 
ment with the results of Fig. |3l 

Finally, we comment on the relation of our scenario 
of domain formation in a an intrinsically nonequilibrium 
system with the standard Kibblc-Zurek mechanism. It 
is clear that some of the key elements of the KZ mech- 
anism, such as the adiabatic-impulse-adiabatic scenario, 
or the "freezing" of correlation length before the phase 
transition, are missing in the polariton condensation sce- 
nario. Nevertheless, in both cases, the scaling laws re- 
sult from the competition of timescales existing in the 
system. Consequently, they will have similar mathemat- 



ical forms. If we assume that the control parameter e 
changes like e ~ {vt) K , the local departure of the or- 
der parameter from the initial "vacuum" state is propor- 
tional to |'i/'| 2 ~ exp(v K t K+1 ), and the quasi-equilibrium 
is achieved at t ~ to a logarithmic factor. 

The correlation length at this point is £ ~ v -k/z{i+k) _ 
These formulas are reminiscent of the standard KZ scal- 
ings ([1]). For example, in the 77 3> 1 limit e(t) ~ t is 
a good approximation over the entire dynamics up to t, 
and we have K = 1 and z = 2 before condensation, which 
results in the correct prediction Nd ~ 7/ 1//4 . 

In conclusion, we investigated the formation of coher- 
ence of exciton-polariton condensate in a one-dimensional 
nanowire. We demostrated the spontaneous forma- 
tion of exciton and polariton domains and found that 
the standard adiabatic-impulse-adiabatic scenario of the 
Kibble-Zurck theory cannot describe the system dynam- 
ics. Nevertheless, power-law critical scalings for the num- 
ber of created defects and the condensation time are 
still present, which results from the competition between 
characteristic timescales in the system. Our study can 
have practical implications for applications proposed for 
polariton systems, as we point to the possible difficulty 
in creating a coherent state on a finite timescale. 
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